arXiv:1509.03449vl [cond-mat.str-el] 11 Sep 2015 


Strong correlations elucidate the electronic structure and phase-diagram of 

LaAlOa/SrTiOs interface 

E. Maniv, M. Ben Shalom, A. Ron, M. Mograbi, A. Palevski, M. Goldstein, and Y. DagarQ 
Raymond and Beverly Sackler School of Physics and Astronomy, Tel-Aviv University, Tel Aviv, 69978, Israel 


The interface between the two band insulators SrTiOa 
and LaAlOa unexpectedly has the properties of a two di¬ 
mensional electron gas. It is even superconducting with a 
transition temperature, Tc, that can be tuned using gate 
bias Vg, which controls the number of electrons added 
or removed from the interface. The gate bias - tempera¬ 
ture (Vg, T) phase diagram is characterized by a dome¬ 
shaped region where superconductivity occurs, i.e., Tc 
has a non-monotonic dependence on Vg, similar to many 
unconventional superconductors. In this communication 
the frequency of the quantum resistance-oscillations ver¬ 
sus inverse magnetic field is reported for various Vg. This 
frequency follows the same nonmonotonic behavior as Tc; 
similar trend is seen in the low field limit of the Hall co¬ 
efficient. We theoretically show that electronic correla¬ 
tions result in a non-monotonic population of the mobile 
band, which can account for the experimental behavior of 
the normal transport properties and the superconducting 
dome. 

INTRODUCTION 

When LaAlOs is epitaxially grown on Ti02-terminated 
{100} SrTiOa, conductivity appears starting from a 
LaAlOa thickness threshold of four unit cells [U [2] . The 
transport properties as well as the superconducting ones 
are strongly dependent on gate bias W- Recently, 
conducting interface has also been observed in {110} 
SrTiOs/LaAlOa interfaces [5]. Gate bias can be either 
applied by biasing the back of the substrate relative to 
the conducting layer (back gate) or by applying the elec¬ 
tric field across the thin dielectric LaAlOa layer (top gate) 
[7] . Superconductivity at the interface is observed over a 
carrier density range of a few 10^^ cm“^ but is suppressed 
at higher densities. In bulk SrTiOs superconductivity 
is known to depend on the number of charge carriers 
[SlIS], and it appears in one or two bands depending on 
carrier density m- In contrast with SrTiOa/LaAlOs, 
in bulk SrTiOa superconductivity extends over a broad 
region of charge densities. It has been proposed that 
two-dimensional fluctuations destroy superconductivity 
in the low carrier regime [3]; However, it is still difficult 
to understand why superconductivity disappears at a rel¬ 
atively low density on the overdoped side of the phase 
diagram of SrTiOa/LaAlOa. 

Previously we have found evidence for the existence 
of multiple bands in SrTiOa/LaAlOa pT] . Santander- 
Syro et al. showed that in conducting SrTiOa surface 


the degeneracy of the titanium (Ti) t 2 g bands is removed 
by the surface, and they split into a lower energy d^y 
band and dy^, dxz bands that are filled at higher gate 
biases [12) . Joshua et al. introduced an atomic spin-orbit 
coupling term that mixes the three bands into a new, 
more complex, band structure |18j . Nernst effect and 
magnetotransport measurements [T3] suggested that the 
mobile band is the one responsible for superconductivity. 
Another view is by Joshua et al. , who suggested that 
a second band becomes populated exactly at a critical 
density where the superconducting critical temperature 
T(. is maximal m- This issue is still a matter of debate. 

In this paper we use the Shubnikov-de Haas (SdH) ef¬ 
fect (probing the area of the Fermi surface), as well as 
the resistivity and the low field limit of the inverse Hall 
coefficient (which is dominated by the density of the 
mobile charge carriers) to study the band structure of the 
SrTiOs/LaAlOs interface. We show that the SdH fre¬ 
quency and follow a nonmonotonic behavior similar 
to the dependence of the superconducting critical tem¬ 
perature Tc and magnetic field He on gate bias. Both 
the SdH frequency and exhibit an anomalous de¬ 
crease upon increasing the gate bias beyond maximum 
Tc. We interpret this unconventional decrease as aris¬ 
ing from electronic correlations between the Ti t 2 g bands 
that are mixed by the atomic spin-orbit interaction. Our 
calculations show that in this case the population of the 
mobile band and its density of states are nonmonotonic 
functions of the chemical potential (/r). This can explain 
the peculiar gate dependence of the transport properties, 
as well as the decrease in Tc on the overdoped regime in 
the (Vg, T) phase diagram. 


RESULTS 
Experimental Data 

We use three samples of {100} LaAlOa/SrTiOa inter¬ 
face patterned into Hall bars (see Methods for full de¬ 
scription). Sample A was used to study the supercon¬ 
ducting dome and the low field Hall effect and Samples 
B and C were used for studying the low field Hall effect 
and the gate dependence of the SdH frequency. 

Figurej^presents the sheet resistance Rn of Sample A, 
normalized by its value at T=800 mK for the same gate 
bias. Defining Tc as the temperature for which R(Tc) = 
R(T=800 mK)/2, the green contour represents Tc as a 
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FIG. 1: Superconducting phase diagram. A color map of 
the normalized resistance plotted as a function of gate bias 
and temperature. Each resistance value is divided by the 
resistance at T=800 mK (above Tc) and the same value of 
Vg (Sample A). 


function of Vg. Tc is maximal for Vg = 2.6 ± 0.3V. The 
entire superconducting dome is revealed by varying Vg 
by merely 15 V, thus demonstrating the performance of 
our gate design m- 

The shape of the superconducting dome is somewhat 
different from the previously reported one [3]. In par¬ 
ticular, in the vicinity of the base of the dome the gate 
dependence of Tc is very sharp. At both edges of the su¬ 
perconducting dome the equi-resistance lines are almost 
perpendicular to the gate bias axes. For Vg < —4.8 V this 
behavior persists into the normal state and hence may 
not reflect intrinsic properties of the material, but may 
arise from the contacts becoming non-ohmic. This region 
is not relevant for our experiment. The fine-tunability of 
our devices allows us to study in detail the dependence of 
the magneto-transport properties and the superconduc¬ 
tivity on gate bias. 

Figure [^a) presents the longitudinal resistivity as 
function of magnetic field for several fixed gate biases 
for Sample A. Superconductivity is suppressed at fields 
greater than the critical magnetic field. He, defined as 
R(Hc)=R(T=800mK)/2. 

In the lower inset of Figure [^b) we show the sheet 
resistance of Sample B versus high magnetic field for 
Vg=3V. SdH oscillations can be seen. In order to make 
them clearer, we subtracted a monotonic polynomial 
background from the data and used a low pass filter to 
remove noise (Figure [^b)). Additional data for Sample 
B is presented in Supplementary Figures 1(b) and 4(a), 
and for Sample C in Supplementary Figures 2 and 4(b). 
The upper inset shows a fast Fourier transform used to 
find the SdH frequency for the data in the lower inset. 
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FIG. 2: Low and High field measurements, (a) Sheet resis¬ 
tance is plotted on a logarithmic scale as a function of mag¬ 
netic field at T=60 mK for various fixed values of gate bias 
Vg (Sample A), (b) Resistance versus inverse magnetic field 
after subtraction of a smooth polynomial background for vari¬ 
ous fixed values of gate bias Vg (Sample B). Successive curves 
are shifted by 517 for clarity. Lower Inset: The resistance per 
square at Vg=3V versus the magnetic field. Upper Inset: Fast 
Fourier Transform (FFT) of the Vg=3V data. 


We found the frequency to be independent of the analysis 
and background subtraction as explained in the Methods 
section. 

We note that Sample B had both back and top gate 
electrodes. The latter was prepared using an additional 
e-beam lithography process, which apparently ruined the 
superconducting transition. However, all other proper¬ 
ties of this sample remained similar to those of samples 
A, C and other samples previously measured. 

We shall now examine carefully the gate dependence 
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of the superconducting and transport properties in or¬ 
der to demonstrate that the low field limit of the Hall 
coefficient, the frequency of the SdH effect, Tc, and He 
all exhibit similar nonmonotonic behavior as a function 
of gate bias. In Figure [^a) the inverse low-field Hall 
coefficient in units of carrier density is plot¬ 

ted as a function of Vg measured at 60mK for Sample 
A. Below IK is practically temperature inde¬ 

pendent. Yellow squares are data taken with a fixed Vg 
while scanning the field from -3T to -I-3T and fitting a 
linear curve up to 3 T to the antisymmetric part. The 
low-field data is completed by the magenta circles, which 
represent the Hall number inferred from the transverse 
voltage measured for /xoH = ±3 T after antisymmetriz¬ 
ing the data. The data taken using these two methods 
coincide, indicating that the Hall signal is approximately 
linear in field up to 3T, as can be seen in Figure [^a). A 
clear non-monotonic signal is observed with a maximum 
at Vg=2.45±0.3 V. We also plot He (blue diamonds), de¬ 
fined by R(He)=R(T=800mK)/2. He, Te and 
are all non-monotonic. The maximum of {eRH)~^ and 
that of Te and He appear at the same gate bias within 
error: Vg=2.6±0.2 V. 

The nonmonotonic behavior of appears to be 

a universal property of the LaAlOa/SrTiOa interface, 
independent of LaAlOs layer thickness (between 6 to 
16 unit cells). It can also be seen in other reports with 
various film thickness: e.g., [1](10 unit cells), [13] (6 and 
10 unit cells), [Il](10 unit cells). In Figurej^a) it is pre¬ 
sented for 6 unit cells with a much higher resolution. This 
behavior is reproduced for 16 unit cells in Figure [^b), 
and for a third sample of 16 unit cells in Supplementary 
Figure 3. 

To better understand which bands govern the low field 
Hall coefficient (ei?//)“^ we studied it together with the 
SdH frequency, which is a direct clean measurement of 
the number of carriers in the mobile band. The SdH 
signal of less mobile bands is strongly suppressed by its 
exponential dependence on the Dingle scattering time. 
We plot (ei?_f/)“^ and the SdH frequency for Sample B 
as function of back gate bias in Figurej^b). Blue squares 
are the Hall coefficients (ei?//)“^ at low magnetic fields 
extracted from a linear fit to the Hall resistance data in 
Supplementary Figure 1(c) between OT and 2T. Green di¬ 
amonds are the SdH frequency, calculated from FFT of 
the data in Figure [^b). Both quantities follow the same 
nonmonotonic behavior, despite an order-of-magnitude 
difference in the carrier densities they correspond to, as 
previously noted by Ben Shalom et al. m]. This means 
that the population of the mobile band decreases with in¬ 
creasing total number of carriers (increasing gate bias). 
The nonlinear response of SrTiOa to electric field m 
cannot account for the nonmonotonic effects that we ob¬ 
serve, since the dielectric function is monotonic as func¬ 
tion of the electric field, and, moreover, the total carrier 
density should in any case increase monotonically with 
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FIG. 3: Non-monotonic behavior of the transport and super¬ 
conducting properties, (a) (Sample A) Left axis (indicated 
by yellow and magenta arrows): The inverse of the Hall co¬ 
efficient is plotted as a function of the gate bias for a fixed 
magnetic field of 3 T (magenta circles). The Hall coefficient 
was also extracted by measuring the Hall resistance as a func¬ 
tion of the magnetic field from 0 T to 3 T and using a lin¬ 
ear fit (yellow squares). Right axis (indicated by a blue ar¬ 
row): Superconducting critical magnetic field He, defined as 
R{T = 60mA, He) = 1/2R{T = 800mK,H = 0), plotted for 
different fixed gates bias values at T=60 mK (blue diamonds), 
(b) (Sample B) Left axis (indicated by a blue arrow): The in¬ 
verse of the Hall coefficient inferred from a linear fit to the 
data in Supplementary Figure 1(c) up to 2T is plotted as a 
function of the gate bias (blue squares). Right axis (indicated 
by a green arrow): The SdH frequency is plotted as a func¬ 
tion of the gate bias (green diamonds). The SdH frequency is 
calculated from FFT analysis of the data in Figure 2(b). 


the applied gate voltage. Therefore, such a scenario is 
possible only if electronic interactions are considered, as 
we now discuss. Sample C showed the same behavior of 
decreasing SdH and (ei?_f/)“^ for increasing gate bias (see 
Supplementary Figure 3). 
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Theoretical Analysis and Results 

To interpret our results we employ a minimal theoret¬ 
ical model of three t 2 g bands (da;j,, and dy^), split 
and mixed by atomic spin-orbit coupling, as suggested 
by Joshua et al. |13| . The main new ingredient in our 
model is a Hubbard-type repulsion between electrons oc¬ 
cupying different orbitals in the same unit cell. We use a 
mean-field approximation in which an electron in a cer¬ 
tain orbital i = xy, xz, yz experiences an average increase 
of UijNj in its energy due to repulsion by electrons in 
orbital j, with Nj being the average occupancy of elec¬ 
trons in the orbital per unit cell (summed over spin). 
The intra-band interaction between electrons of opposite 
spins is \UiiNi. The orbital populations are calculated 
self consistently, as explained in the theoretical methods 
section. For simplicity we took the interaction energy 
Uij = U = 2.7eV both for the intra- and inter-band in¬ 
teractions for all orbitals, consistent with Ref. [18]. 

Due to the interactions, upon increasing the chemi¬ 
cal potential (i.e., increasing the gate bias) the band 
structure is modified, as demonstrated in Figure Qa)- 
(c), where the band structure is presented for different 
values of y. Each of the three bands should be slightly 
split due to the Rashba spin-orbit coupling. We neglect 
this effect since it is very small, as demonstrated in Sup¬ 
plementary Figure 6 (although it is important for the 
parallel-field dependence of superconductivity [5| and for 
the weak anti-localization effect [19|). 

In addition, scanning SQUID measurements detected 
free magnetic moments at the interface |20j and magne¬ 
tization measurement [211122j found evidence for a mag¬ 
netic phase at elevated temperatures. However, the lat¬ 
ter measurements were performed at significant magnetic 
fields, presumably sufficient to align the moments. We 
have recently found evidence for a magnetic order at zero 
field below IK [23]. The energy scale of this interaction 
is very low compared with the terms in the presented 
Hamiltonian. Therefore, and in order to keep our calcu¬ 
lations simple, we did not allow for a Stoner instability. 

The resulting total densities of electrons in each band 
are presented in Figure Qd). While the electron density 
in the lowest-energy band 1 [red line in Figure|^a)-(c)] is 
monotonically increasing with /i, the electron density in 
band 2 which is populated next [blue line in Figure |^a)- 
(c)] has a nonmonotonic dependence on /r (note that the 
overall electron density is still monotonic with y). In 
Figure l^e) we present the calculated single-particle den¬ 
sity of states (DOS) for band 2, which exhibits a similar 
nonmonotonic /i-dependence. 

To calculate the Hall resistance we employ a model of 
point impurities. In the absence of orbital mixing this 
would lead to momentum-independent scattering rates 
for each band. In that case, in the two-band regime, 
the low-field limit of the inverse Hall coefficient (ei?i/)“^ 


would be (eRn) ^ , with Q being the mo- 

bility and the carrier density of the unmixed Hh band 
(orbital). With band mixing the situation is more com¬ 
plicated, since the composition of the bands in terms 
of the original t 2 g orbitals changes as function of the 
momentum and as function of /i. We performed the 
corresponding calculations, taking the carrier densities 
from Figure l^d), and assigning mobilities of 6000 and 
500 cm^/V-s to the unmixed d^y and dyzjdxz bands, re¬ 
spectively, roughly representing the ratio of the unmixed 
bands’ effective masses. The resulting low-field inverse 
Hall coefficient (ei?ij)“^ is shown in Figure ge). The 
calculated {eRn) ^ is nonmonotonic as function of the 
chemical potential, and its maximum appears at the same 
H as the maximum of the DOS of band 2 (blue). This 
behavior is similar to the observed dependence of the 
measured (ei?//)“^ as function of gate bias for all sam¬ 
ples. 

DISCUSSION 

In a simple capacitor one expects the carrier concen¬ 
tration to monotonically increase with gate voltage (even 
in the presence of nonlinear dielectric materials such as 
SrTiOs). By contrast, in our experiment we see a non¬ 
monotonic behavior of the SdH frequency and (ei?jj)“^ 
as function of gate bias. Our theoretical model does give 
rise to such nonmonotonic behavior as seen in Figure [^ 
What is the reason for this theoretical result? Electronic 
interactions cause a competition between the occupan¬ 
cies of different bands: When the population of one band 
increases, the energies of the states in the other band are 
shifted and the other band’s occupation decreases. In 
the presence of such a competition occupying the band 
with the highest single-particle density of states is usu¬ 
ally preferable, since this results in lower single-particle 
energy. As Figure Qa)-(c) reveals, at low energies band 
1 is mainly composed of low effective mass (low density 
of states) dxy orbitals, but at higher energies it gains a 
heavier (high density of states) dxz/dyz character; the 
opposite occurs for the band 2. Thus, as the gate bias, 
and hence /r, are increased from large negative values, the 
following occurs: (i) Initially band 1 is occupied, since its 
energy is lower. As it becomes progressively filled, the en¬ 
ergy of band 2 increases [Figurej^a)]; (ii) At some point 
/r reaches the bottom of the shifted band 2, where band 
2’s density of states is large, and it begins to fill up [Fig- 
urej^b)]; (hi) At higher /i, when band 1 becomes heavier 
than band 2, it becomes preferable to increase the pop¬ 
ulation of band 1 further at the expense of band 2 [Fig¬ 
ure l^c)] . The population of band 2 is thus nonmontonic 
as function of y. A similar population switching effect 
has recently been discussed in the context of quantum 
dot physics, where a narrow (high density of states) level 
may increase in occupancy at the expense of a broad (low 
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FIG. 4: Results of the theoretical model, (a-c) Band structure including the effects of mean-field electronic correlations for 
three different chemical potentials /r. A cut through the center of the Brillouin zone along the kx axis is shown. The calculated 
carrier densities for each band are shown, (d) The carrier concentration of the two lower bands (band 1 and band 2). A 
non-monotonic behavior of the carrier density of band 2 is obtained, (e) Left axis (indicated by a yellow arrow): Calculated 
inverse Hall coefficient in the low field limit as a function of /r. Right axis (indicated by a blue arrow): The density of states 
per unit cell of band 2. Both properties show nonmonotonic behavior with a maximum at a similar value of fi. 


density of states) one (See, e.g.. pi^B5] ). 

In Figure]^ a) we compare the calculated carrier den¬ 
sity of band 2 (blue) and the one extracted from the 
frequency F of the SdH signal for Sample B through 
n = eF/irh. First, we note that the calculation gives 
the right order of magnitude for the carrier density. Sec¬ 
ond, the nonmonotonic behavior of the SdH frequency 
is reproduced by the calculation. The chemical potential 
was shifted and scaled such that the maximum in the cal¬ 
culated carrier density and the maximal measured SdH 
frequency coincide. This scaling is consistent with the 
estimated capacitance of our devices m- Furthermore, 
Since both bands contribute to the Hall coefficient, while 
only the mobile second band (whose population is sup¬ 


pressed by the interaction effects) contributes to the SdH, 
our model naturally explains the order of magnitude dif¬ 
ference between (ei?ij)“^ and the density corresponding 
to the SdH frequency, which has hitherto been a puzzle 

m- 

Our analysis implies that band 1 has much lower mobil¬ 
ity than band 2. This could be due to the fact that in the 
gate voltage region corresponding to the nonmonotonic 
behaviors band 1 is mainly composed of the heavier dxz- 
dyz orbitals, whereas band 2 mainly includes the lighter 
Axy orbital. It should be noted that at lower values of /r 
(before band 2 enters) band 1 has a dxy character, but is 
still not observed in SdH. This can be accounted for using 
the density functional theory calculations of Belugas et 
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FIG. 5: Comparing experimental data and theoretical calcula¬ 
tions. (a) Bottom-Left axes (indicated by a blue arrow): The 
carrier density of band 2 [taken from Figure 4(d)] is plotted as 
a function of /i. Top-Right axes (indicated by a green arrow): 
The carrier density inferred from the SdH frequency in Figure 
3(b) is plotted as a function of the gate bias Vg. Both show 
a similar behavior as function of Vg or fj, (which are related 
through the gate capacitance), (b) Top axis (indicated by a 
purple arrow): Critical temperature of Sample A (taken from 
Figure 1) is plotted as a function of Vg. Bottom axis (indi¬ 
cated by a blue arrow): Critical temperature calculated from 
the weak-coupling BCS formula using the density of states 
of band 2 [taken from Figure 4(e)] is plotted as a function 
of chemical potential fi. Both critical temperatures exhibit a 
nonmonotonic behavior with a maximum at the same jj, (Vg) 
as the inverse Hall coefficient. 


al. [30]. They find that at low values of /i the dxy orbital 
resides in Ti atoms at the STO/LAO interface, where 
scattering centers are expected to reside, while at higher 
values of /r the d^y orbital moves to Ti atoms away from 
the interface. This picture is also in line with the trans¬ 
port measurements of Bell et al. |3| , and with the observed 


difference between the effects of top and bottom gates on 
the transport properties [7|. These two gate configura¬ 
tions can effect the location of the wave function along 
,—, the z direction in an opposite way. 

'B Since the theoretically-calculated DOS p 2 of the band 2 

'—' (blue) has a similar nonmonotonic behavior as the mea- 

sured (ei?//)”^ and Tc, it is tempting to try and com- 
S pare the measured Tc with the weak-coupling single-band 
Q BCS formula, Tc « In order to obtain 

.H the correct order of magnitude for Tc with the measured 

§ Debye temperature Td = dOOK of SrTiOa [3T], we use 
electron-phonon coupling energy Vbcs = 0.196 eV. This 
value is within the weak coupling limit for all the calcu¬ 
lated DOS values used. 

In Figure [^b) we plot the calculated Tc together with 
the measured one (taken from Figure with /i and 
Vg scaled to match the superconducting dome region. 
The nonmonotinc behavior of Tc is reproduced. Using 
the same scaling the calculated and measured (ei?//)“^ 
showed the same behavior as a function of p, and Vg (see 
Supplementary Figure 5). While using a single-band BCS 
model is probably an oversimplification (and indeed for 
the low carrier densities a non-BCS type behavior has 
been reported in planar tunneling measurements [32]), it 
is interesting to see that it can nicely capture the behav¬ 
ior of Tc versus gate bias in a SrTiOa/LaAlOa interface, 
and in particular account for the decrease in Tc on the 
overdoped regime. 

In Summary, we found that the Shubnikov-de Haas 
frequency, the inverse low field Hall coefficient and the 
superconducting transition temperature all exhibit a sim¬ 
ilar nonmonotonic dependence on gate bias. In order to 
explain our data we employed a model involving elec¬ 
tronic correlation between the spin-orbit split titanium 
t 2 g bands. In this model the second band, which is pop¬ 
ulated in a non-monotonic fashion at higher gate bias 
values, is the one responsible for the SdH effect and for 
superconductivity. 

METHODS 
Experimental 

We grow epitaxial layers of LaAlOa using reflec¬ 
tion high energy electron diffraction (RHEED) monitored 
pulsed laser deposition on atomically flat Ti 02 termi¬ 
nated {100} SrTiOa 0.5mm thick substrates in stan¬ 
dard conditions, oxygen partial pressure of 1 • 10“^ Torr 
and temperature of 780°C', as described in [33] ■ In the 
first step two unit cells of LaAlOa are deposited. Then 
Hall bars of 9 x 3 pm? (A) 12 x 3 pm^ (B) and 20 x 5 
pm^ (C) are patterned using electron beam (A) or op¬ 
tical (B and C) lithography followed by deposition of a 
40 nm thick amorphous BaTiOa layer and lift-off to de¬ 
fine the conducting channel, where an additional layer of 
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LaAlOa is epitaxially grown in the final step. The total 
LaAlOa thickness is 6 (16) unit cells for Sample A (B 
and C). The design aims to minimize screening from the 
contact pads and leads, thus enhancing the gate bias ef- 
fectivity. The strong response to back-gating is possible 
for dielectric substrates with high permittivity (high e) 
when the distance d between the gate electrode and the 
tunable surface (or interface) becomes much larger than 
the width of the mesoscopic conducting channel i. This 
results in a capacitance per area of the order of 5 x 10^^ 
[electronsxcm“^Volt“^] for d = 0.5 mm, £=3 /rm, an 
order of magnitude larger than for standard planar ca¬ 
pacitor geometry m- 

Gold gate electrodes are evaporated to cover the back 
of the substrate. The leakage current is unmeasurably 
small (<lpA). The resistivity and low field measurements 
of Sample A were performed in a dilution refrigerator 
with a base temperature of 60 mK using a Princeton Ap¬ 
plied Research 124A lock-in amplifier. During cool down, 
the gate bias was scanned back and forth while recording 
the temperature and measuring the sheet (per-square) 
resistance. The gate scan is reversible within a 2 mV 
resolution, as long as the maximal gate bias (10 V) is not 
exceeded. The high field measurements were performed 
on Sample B and C. Sample B was measured in a di¬ 
lution refrigerator with a base temperature of 20mK at 
magnetic fields of up to 18 Tesla and a constant top gate 


voltage of 1 V. Sample C was measured in a He3 cryostat 
with a base temperature of 350mK and at magnetic helds 
of up to 34.5 T. Both samples B and C were measured 
using a Lakeshore 370 resistivity bridge with 3716L low 
resistance preamplifier and scanner at Tallahassee Na¬ 
tional High Magnetic Field Laboratory. 

For the SdH analysis we first subtracted a smooth poly¬ 
nomial background as described in ref m- We made 
sure that the background consists of no oscillations and 
that the frequency is independent of background chosen. 
The background was fitted using a polynomial with a 
maximum order of 4. The frequency is found both by us¬ 
ing FFT analysis with a step-type window function and 
manually (using the distance between the extrema in the 
data). Both methods give similar results within an error 
margin that we estimate to be 5-10%. 

Theoretical 

We describe our system using the model introduced 
in Ref. [T3], to which we add local (Hubbard) interac¬ 
tions in the mean-field approximation. In k-space the 
Hamiltonian is a 6 x 6 matrix in the basis of 3 orbitals 
(i = xy,xz,yz) 0 2 spin states (s =t,4,). It is a sum of 
three terms, H{k) = i7o(k) + -ffso(k) + i7int(k), the first 
of which being 


with 


i?o(k) 


( 0 exz(k) 

y 0 2 td sin(fc 2 ,a) sin(fcj,a) 


0 


2td s\Ti{kxa) sin(fcya) 


^liz (k) 


® I 2 


ea;y(k) =2t/ [2 — cos{kxa) — cos(A:ya)] 

£ 3 , 2 (k) =2ti [1 - cos{kxa)] + 2th [1 - cos(fcya)] 
£y^(k) =2th [1 - cos(fc 2 ;a)] -I- 2ti [I - cos(fcya)], 


( 1 ) 


( 2 ) 

(3) 

(4) 


where a = 3.905A is the SrTiOs lattice constant, ti = 
875meV, th = 40meV, and td = 40meV are, respectively, 
the light, heavy, and yz-xz mixing tunneling matrix ele¬ 
ments. Through the relation ti/h = ^^/(2m;/;ia^), these 
correspond to light/heavy effective masses of mi = 0.7me 
and mh = 15me, respectively. Ae = 47meV is the rela¬ 
tive shift between the xy and xzjyz orbitals due to the 
confinement potential at the SrTiOa/LaAlOa interface. 
We neglect orbital coupling due to the asymmetric con¬ 
fining potential and the Rashba term (see Supplementary 
Figure 6 ). Following XAS [31], ARPES [T3], and trans¬ 
port m results, we assume that the other dxy sub-band 
predicted by DFT calculations |30| are localized and do 


not play a role in transport, following in this respect pre¬ 
vious theoretical studies naiis]. 

Atomic spin orbit coupling contributes 

■A ( ^ ^y\ 

ggo(k) = l-CTx 0 CTz j , (5) 

V cTy -a, 0 / 

with Aso = dOmeV. 

Finally, the mean-field interaction term is 

0 0 \ 

Hi„t= 0 E-ft 0 0/2, (6) 

V 0 0 EZj 





where /2)Nj {i,j = xy,xz,yz), 

where Nj is the average occupancy of orbital j (summed 
over spin) per unit cell. The subtracted term in E^^^ re¬ 
moves the unphysical self-interaction of electrons in the 
same orbital and spin state. The matrix [/ij is sym¬ 
metric. Moreover, the symmetries of our system dic¬ 
tate that Uxy^xz — Uxy^yz and Uxz,xz — Uyz^yz- fa 
our calculations we take for simplicity all the element 
to be equal, Uij = U = 2.7eV. The average occupan¬ 
cies Nj are determined self-consistently, taking into ac¬ 
count the contributions of the different eigenstates of the 
Hamiltonian weighted by their decomposition in terms 
of the orbitals and the Fermi-Dirac distribution, /(e) = 
-I- 1]: Denoting the eigenvectors of the ma¬ 
trix H{k) by i/)m,r(k; i, s), with corresponding eigenvalues 
£m,T(k) {i = xy,xz,yz and s denote the bare or¬ 

bital and spin, whereas m = 1, 2, 3 denote the bands and 
r = ±1 is a Kramers spin index), we have 


E 


’^dk.2 


IBZ 


( 2^)2 


IV'm,r(k;j,s)| /[em,r(k)], (7) 


where the integration is over the first Brillouin zone. 

After the self-consistency equations have been solved, 
the total density of electrons populating band m = 1, 2, 3 
(related to the frequency of the SdH oscillations) can be 
calculated as: 





dk 


f [em,r(k)] . 


( 8 ) 


The local density of states (per unit cell) Pm of band m 
at the Fermi energy is the derivative of a^rim w.r.t. p, 
keeping the mean-field interaction-induced energy shifts 
El^^ constant. 


* yodagan@post.tau.ac.il 
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